Correlated electrons systems on the Apollonian network 



O 
O 

(N 

^— > 
O 

O 



3 



C 

o 
o 



O 

-i— > 

a 

E 

T3 ■ 

c : 
o . 
o ■ 



X 



Andre M. C. Souza 1,2 and Hans Herrmann 3 ' 4 
1 Institut fur Computerphysik, Universittit Stuttgart, Pfaffenwaldring 27, 70569 Stuttgart, Germany 
2 Departamento de Fisica, Universidade Federal de Sergipe, 49100-000 Sao Cristovao-SE, Brazil 
3 Institut fur Baustoffe, ETH Honggerberg, HIF E 12, CH-8093 Zurich, Switzerland and 
4 Departamento de Fisica, Universidade Federal do Ceard, 60451-970 Fortaleza-CE, Brazil 

(Dated: February 6, 2008) 

Strongly correlated electrons on an Apollonian network are studied using the Hubbard model. 
Ground-state and thermodynamic properties, including specific heat, magnetic susceptibility, spin- 
spin correlation function, double occupancy and one-electron transfer, are evaluated applying direct 
diagonalization and quantum Monte Carlo. The results support several types of magnetic behavior. 
In the strong-coupling limit, the quantum anisotropic spin i Heisenberg model is used and the 
phase diagram is discussed using the renormalization group method. For ferromagnetic coupling, 
we always observe the existence of long-range order. For antiferromagnetic coupling, we find a 
paramagnetic phase for all finite temperatures. 
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I. INTRODUCTION 

Scale-invariant networks have been subject of inten- 
sive study in view of possible insights into inhomo- 
geneous problems such as random magnets, surfaces, 
porous rocks, aerogels, sponges, etc. They are relevant 
to many different real complex situations like biological, 
social and technological systems 0, Q , where the struc- 
tures can embody a particular graph, e.g., a scale free 
network. In this case, the fraction of sites with k con- 
nections follows a power law [jj . This topology has been 
found, in particular, on World Wide Web and Internet 
networks |4(. 

More recently, the area of scale-invariant networks 
has been highly motivated by the creation of a syn- 
thetic nanometer-scale Sierpinski hexagonal gasket, a 
self-similar fractal macromolecule Q . New types of pho- 
toelectric cells, molecular batteries and energy storage 
may be possible. This perspective is relevant in view of 
the renewed interest in correlated electron systems on 
these networks, like studies of quantum magnetism, su- 
perconductivity, metal-insulator transition, etc. An im- 
portant aspect that can be analyzed is the influence of 
the topology and the appearance of magnetic order. 

In this paper, we investigate one family of free scale- 
invariant graphs, the Apollonian network [ESQ- The 
purpose is to examine correlated electrons on Apollonian 
networks. We explore the connection between the scale- 
invariant network topology and the magnetic properties 
of quantum magnetism models, namely, Hubbard and 
anisotropic spin h Heisenberg models. We have exam- 
ined the Hubbard model using the small cluster diagonal- 
ization 

HQS and quantum Monte Carlo jUGS- The 
Heisenberg model is studied within a real-space renor- 
malization group framework |l3l Il4. Il5j | . 

The construction of the Apollonian network has the 
topology of the contacts of an Apollonian packing of two- 
dimensional disks. As initial configuration we use, for an 
Apollonian packing, three mutually touching circles in- 



scribed inside a circular space. The interstices of the ini- 
tial disks define a curvilinear triangle to be filled. In the 
next iteration, four disks are inscribed, each touching all 
the sides of the corresponding curvilinear triangle. This 
process is repeated indefinitely by setting disks in the 
newly generated curvilinear triangles. The Apollonian 
packing is a fractal whose dimension has been estimated 
as 1.3057 H. 




FIG. 1: First three generations of the Apollonian network. 
Sites represented by squares, diamonds and circles are intro- 
duced in the first, second and third generation. 

The Apollonian network (AN) is derived from the two- 
dimensional Apollonian packing associating vertices with 
the circles and connecting two vertices if the correspond- 
ing circles touch. Fig 1 shows the first three generations 
of the AN. With this construction procedure one obtains 
a deterministic scale free network where the number of 
sites at iteration n is (3" _1 +5)/2. Besides being a scale- 
free network, the AN has interesting properties like being 
Euclidean, matching and space-filling 6J. 

The organization of this paper is as follows. The 
models are briefly introduced in section II. The results 
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for ground-state and thermodynamic properties of the 
Hubbard model on AN are presented in section III and 
IV, respectively. The phase diagram of the quantum 
anisotropic Heisenberg model is discussed in section V 
and conclusions are presented in section VI. 



II. MODELS 

Our study of the magnetic properties on AN are based 
on the Hubbard model. The Hamiltonian of the Hubbard 
model is defined by 



Ti-u = —t 



where 



E ■ 

<ij>a 



Cj a 



(1) 



c'i a i Ci a and rii a = c ia Ci a are respectively the cre- 
ation, anihilation and number operators for an electron 
with spin a in an orbital localized at site i\ the < ij > 
sum runs over nearest-neighbor sites on AN |16|| . 

The question of magnetic order in the one-band Hub- 
bard model has been investigated by several authors and 
much controversy has arisen. In the strong-coupling limit 
a major part of the ferromagnetic phase is predicted. In 
this case, at half filling band, U is much larger than t, 
and the Hubbard model, using a suitable expansion in 
perturbation theory, is formally equivalent to the antifer- 
romagnetic Heisenberg model. The Heisenberg exchange 
parameter J is written in terms of the Hubbard model 
parameters as J = —4t 2 /U, where the Hamiltonian of 
the Heisenberg model is defined by 



He 



-J 



E 

<ij> 



SiSj 



(2) 



where d, is the total spin | operator for the ith site. 
Here we study the anisotropic Heisenberg model to com- 
pare the behavior for different values of the anisotropy 
parameter [l5j |. 
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TABLE I: Spin of the ground state as function of occupation 
number and U/t on the AN of 7 sites. 
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FIG. 2: Ground state spin-spin correlation function versus 
U/t for the half-filled band on the AN of 7 sites. 



III. GROUND-STATE PROPERTIES 

We have obtained exact numerical results of the Hub- 
bard model defined on second generation of the AN, cor- 
responding to seven sites. We computed all the eigenval- 
ues and eigenvectors of the Hamiltonian of Eq. 1 on a 
basis of states for which the occupation number is diago- 
nal. We consider subspaces of fixed total azimuthal spin 
operator Sz- Thus, the maximum possible dimension of 
the matrix to diagonalize is 1225. 

Table I shows the results for the spin of the ground 
state. Until the half-filled band, the ocurrcncc of a ferro- 
magnetic state, where the total spin S is not minimum, 
can be found only for four electrons. In this case the 
ground state is threefold-degenerate corresponding to a 
triplet state of S = 1. This behavior can be explained 
considering the spectrum of the free electron system. Ta- 
ble II lists the single electron energies and degeneracies. 



The double degeneracy in second level forms a triplet in 
lowest energy for the four electrons case. For two, three, 
five, six and seven electrons the ground state always has 
minimum spin. We observe that the Coulombian interac- 
tion does not favour the occurrence of a high spin ground 
state. If U/t is strong, the jumps of electrons decrease 
and a frustrated ordered antiferromagnet is favored. We 
easily see that the network is not bipartite, so its struc- 
ture is antiferomagnetically frustrated. We find a similar 
competition, between interaction and frustration, in the 
low-temperature antiferromagnetic state on a triangular 
lattice [l7| . 

Above the half-filled band we obtain a ferromagnetic 
ground-state. Nagaoka demonstrated that the ferromag- 
netism is expected for the antiferromagnetically frus- 
trated structures in the half-filled band case with one 
excess electron and U — > oo 18]. Here, this limit repre- 



3 



Ei 


Ei — E3 


Ei 


E5 — Ea 


E 7 
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-0.6180 
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1.6180 


1.7525 



TABLE II: Energy levels E t (i=l, 2,..., 7) of free electron 
system (U/t — 0) on the AN of 7 sites. 

sents the eight electrons case for U/t >> 1. We can see 
that if U/t increases the alignment of spins also increases, 
in such a manner that, for U/t > 20.1, the ferromagnetic 
state with the maximum total spin is the ground state. 
For more than eight electrons ferrmagnetism is possible, 
but not stronger than in the case of completely saturated 
ferromagnetism. 

In general, the behavior of short (or long)-range order- 
ing can be better observed studying the spin-spin corre- 
lation function. Considering the azimuthal spin operator 
on site i as Sf = — n^, the spin-spin correlation func- 
tion is defined, for N sites, as 

L s = ^ S*Sf +s . (3) 

i 

The quantity Lq, called local moment, depends on the 
magnitude of the difference between up and down elec- 
tron spin at each site. It shows the degree of localization 
of electrons. For a completely localized system in which 
each site is occupied by a single electron (up or down), 
Lq = 1/4; while for a non-interacting system Lq = 1/8. 
Lg (5 ^ 0) is the correlation between the electron spins 
at different sites. They are related to the magnetic or- 
dering. The results for the ground state spin-spin corre- 
lation function versus U/t in the half-filled band case are 
presented in Fig. 2. L gradually increases if U/t > in- 
creases, i.e., electrons gradually localize. L\ has negative 
sign and increases if U/ 1 > increases and Li has a crit- 
ical value for U/t = 8 where the behavior changes from 
negative to a positive sign. For U/t > 8, L2 is positive, 
L\ is negative, and an antiferromagnetic order appear. 
For U/t = 11.1 there is a change in the spatial symmetry 
of the ground state. It is important to observe the discon- 
tinuous change of L\ and Li. This quantum transition 
point at U/t =11.1 must not be a Mott metal-insulator 
transition because Lq is continuous 0. For U/t >> 1, Lq 
is approximately equal to 1/4 and we can conclude that 
each site is occupied by just one electron. Our results in- 
dicate that the system can present a fcrrimagnetic order, 
however, because the size of the network considered, the 
fcrrimagnetic state has minimum spin S = 1/2, and we 
cannot have high spin ground state. For U/t < 8, L 2 is 
negative and the magnetic behavior of the ground state 
is more complex. 

We analyse the wave function of the ground state, 
defining Pk = | SI/ a,o | 2 , where |^fco| is the nth compo- 
nent of the ground eigenstate. We define the total an- 
tiferromagnetic configuration (k = AF) such that sites 
introduced in the first generation on AN (represented by 
squares in Fig. 1) are occupied by electrons of a type 
of spin and sites introduced in second generation (rep- 
resented by diamonds in Fig. 1) are occupied by elec- 



trons of opposite spin. We observe that Paf/Pu > 2 for 
U/t = and increases to Paf /Pk > 15 if U/t = 8, for all 
k ^ AF, i.e., the total antiferromagnetic configuration 
Paf is approximately 2 times for U/t = and 15 times 
for U/t = 8, more probable than any another configura- 
tion. Thus, we can construct a physical picture where the 
ground state has an antiferromagnetic order (frustrated) 
at any U /t. 



IV. THERMODYNAMIC PROPERTIES 

On one hand, we study the thermodynamic properties 
of the half-filled band Hubbard model using the method 
of small-cluster exact-diagonalization calculations in the 
grand canonical ensemble [ljj. We have calculated all 
eigenvalues and eigenfunctions for AN of 7 sites. On the 
other hand, we study lattices of 63 and 124 sites usin g th e 
grand canonical quantum Monte Carlo method [Til Il2|. 
We used a discrete Hubbard-Stratonovich transformation 
to convert the problem into one of free particles interact- 
ing with a time-dependent Ising field, together with an 
exact updating algorithm for the fermions Green's func- 
tion to compute the relative weights of the Ising config- 
urations. 

We consider several values of U /t for all the stud- 
ied thermodynamic properties. We determine with con- 
siderable accuracy the temperature dependence of the 
spin-spin correlation functions, spin susceptibility, spe- 
cific heat, double occupancy and one electron transfer. 

The temperature dependence of Lq, Li and Li for some 
typical values of U/t is shown in Figs. 3a-3c. If U/t 
increases, Lq gradually increases for U/t > 0, indicat- 
ing that electrons are gradually localizing. At high tem- 
perature Lq gradually decreases and electrons gradually 
delocalize. At very low temperatures, the temperature 
dependence with negative signs of L\ induces an antifcr- 
romagetic ordering. The inset in Fig. 3b shows Monte 
Carlo results for AN with 7, 43 and 124 sites and U/t = 1. 
Results for 43 and 124 sites do not reveal any new be- 
havior, however, the antiferromagetic ordering increases. 
We observe that the properties obtained for 7 sites must 
be qualitatively equivalent to those for larger AN. We ob- 
serve a competition between interaction and frustration 
in the ordered antiferromagnetic state by the analysis of 
Li . For small or intermediary U/t, the negative sign of 
Li is a result of the frustrated structure of AN, disfavor- 
ing the antiferromagnetic order. For U/t >> 1, Li has 
positive sign, showing that the strong electron interaction 
favors the antiferromagnetic order. 

Fig. 4 shows the temperature dependence of the 
magnetic susceptibility, \- The Curie- Weiss behavior, 
X oc 1/(T — 9), is observed for high temperatures. For 
small and intermediary U /t the curves present a peak. 
The temperature associated with the peak of x corre- 
sponds approximately to the rapid decay of L\ and Li, 
and is a consequence of the collective excitations that 
lead to the destruction of the antiferromagetic order [j| . 
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FIG. 4: Magnetic susceptibility as a function of temperature 
with U/t = 0, 0.5, 4, and 8 for the AN of 7 sites. Inset: 
Magnetic susceptibility with U/t — 1 for 7, 43 and 124 sites 
on AN using Monte Carlo method. 
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FIG. 3: Correlation functions Ls vs. temperature with U/t = 
0.1, 1, 6 and 8 for the AN of 7 sites, (a) L ; (b) Li; (c) L 2 . 
In the inset of (b), L\ with U/t — 1 for 7, 43 and 124 sites on 
AN using Monte Carlo method. 



Here the ground state has spin S = 1/2 and x must go 
to infinity at temperature T = 0. The susceptibility x 
with U/t = 1 for 7, 43 and 124 sites on AN using Monte 
Carlo method is shown in the inset of Fig. 4. In this 
case, we do not see any discrepancy between the results 
for different clusters. 

The temperature dependence of the specific heat is 
shown in Fig. 5 for typical values of U /t. For small 
values of U/t there is a peak in the specific heat. In- 
creasing U/t, the peak splits into two, which reflects a 
rearrangement of the fermionic structure in the system. 
The low-temperature peak arises due to low-lying collec- 
tive excitations, while the high-temperature broad peak 
comes from single-particle excitations. This behavior is 



FIG. 5: Specific heat C/Nks as a function of temperature 
with U/t = 0,0.5, 4 and 8 for the AN of 7 sites. 



qu ite g eneral and has been noticed for different structures 
9, 10]. 

Next, we study the thermal average of the double occu- 
pancy d = -k n il n il an d °f the- one-electron transfer 
defined by 

j - E 4 a c ja = U -^L. (4) 



<ij>a 



t 



Fig. 6 shows the temperature dependence of j and d 
(inset) with U/t = 0.1, 1 and 4 for the AN with 124 
sites. Clearly these functions are related to the local 
moment. It is easy to show that for the half-filled band 
d = (1 — 4L )/2. If U/t increases, the local moment 
increases (see Fig. 3a) and the double occupation and 
one-electron transfer decrease. In the strong coupling 



5 



■O 1,2 




FIG. 6: One-electron transfer j as a function of temperature 
with U/t = 0.1, 1 and 4 for the AN of 124 sites. Inset: Double 
occupancy d as a function of temperature. 




FIG. 7: Renormalization group transformation associated to 
the AN. 



limit (U/t >> 1) each site is occupied by just one electron 
and the electronic itinerancy vanishes. This region is 
characterized by localized magnetism and we can apply 
the Heisenberg model. 



V. LOCALIZED MAGNETISM: 
RENORMALIZATION GROUP APPROACH 

The fractal and self-similar nature of scale-invariant 
networks, generated through a recursive procedure, offer 
an appropriated way to get rigorous result at the phase 
transition. The difficulty with the spin-decimation pro- 
cess which leads to the renormalization-group equation 
as found on Bravais lattices due to parameter prolifer- 
ation, disappears on scale-invariant networks pel l'2ll| . 
Due to the fact that fractal lattices are characterized by 
dilation invariance but do not have translational invari- 
ance, the study on these networks can with restriction be 
used to model Bravais lattices. Particularly self-similar 
networks can mimick Bravais lattices for some magnetic 
models, providing equal critical temperatures and expo- 
nents [2ltl2^ . 

The real-space renormalization group (RG) approach 
has been applied with success to the study of the 



anisotropic spin t, Heisenberg model on self-dual hier- 
archical lattices pLia. fl4L Il5j| . Here, we apply this analysis 
on AN. It is defined by the cluster transformation pro- 
cess, preserving the Hamiltonian form, as illustrated in 
Fig. 7. 

The dimensionless Hamiltonian is defined by 
4J 



-PH 



<ij> 



[(l-A)(SfS? + S?S1) + S!Sfi,{5) 



where (3 = l/kgT, (ij) denotes first-neighboring lattice 
sites, A is the anisotropy parameter and the Sf{a = 
x,y,z} is the cv-spin 1/2 operator on site i. The RG 
recurrence equation is obtained by imposing 



exp(/i 123 + C) 



Tr exp(fti 234 ) 

site 4 



(0) 



where A1234 and h 123 are, respectively, the Hamiltonian 
of the four-site cluster and of the renormalizcd three-site 
cluster, shown in Fig. 7. To make the RG equation pos- 
sible, an additive constant C has been included. The 
RG Eq. 6 establishes the relation between the set of pa- 
rameters (J, A) and the set of renormalized parameters 
(J , A ). The non-commutativity between the Hamilto- 
nians associated with neighboring clusters is neglect, and 
therefore our results are approximations for all tempera- 
tures, only being asymptotically exact at high tempera- 
tures PH. 

Defining K = J/ksT and expanding exp(h 123 ) as 



+S2SI 
we obtain that 



+ 4b (S^SZ + SfS% 
SZSi)+4c(SlS z 2 - 
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c 

C = -3k 
Analogously, 
exp(/ii234)= a 
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1 c -a 
T 

o c — a 



2b 
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a 



2b'' 



+ ln(a -I- 3c ). 



(8) 



> + 4 £ (i<i) MSfS* + S?S?) + cSfSI] 

+16 J2(i <3 )=£(k<i) [d(SfSj + SfSj)S%.S( 

+g(SfS* + S?S*)(S* k S? + SIS?)] 

+16 / S{S^Sl (9) 

where a,b,c,d,g and / are functions of K and A. It 
is easy to see from Eq. 6 that a —2 a, b = 26 and 
c = 2 c. The set of parameters (K , A ) as functions 
of (K, A) can be determined by diagonalizing h\ 2 z and 
/11234 and using Eqs. 6-9. After some calculus we obtain 
the analytical RG equation as 
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FIG. 8: Flow diagram for the Apollonian cell of Fig. 7. Open 
circles, and full squares respectively, denote, the semi-stable 
and fully stable fixed points, (a) ferromagnetic case (J > 0). 
(b) antiferromagnetic (J < 0) case. 
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where K , A are functions of c a , b given by Eq. 8. 
We observe that the Ising and isotropic Heisenberg mod- 
els are mapped into themselves. The RG recurrence is 
simplified as 



1 1 + e 6K 

k = — ln(- tttt) 

4 K l + e~ 2K ' 



and 

k = 



1 3 + 5e SK 
6 n ^6 + 2e- 4K ' 



(11) 



(12) 



corresponding to the Ising (A = 1) and the isotropic 
Heisenberg (A = 0) limits, respectively. 

Fig. 8 shows the phase diagram in the (fc^T/ J, A) 
space. In the ferromagnetic case (fcsT/ J>0), we always 
observe the existence of ferromagnetism(F) independent 
of A. This result has been found on AN for the Ising 
model using the transfer matrix method and on other 
scale- free lattices \lQ . Those results agree very well with 
ours. Here, we observe that the quantum fluctuations of 
the XY part of our Hamiltonian does not destroy this 
ordering generated by the topology of the AN. In the 
antiferromagnetic case (fc^T/J < 0), we verify the non- 
existence of the ordered phase for all A at finite T. Again, 
the results obtained are similar and consistent with the 
Ising limit 0|. 



VI. CONCLUSIONS 

In conclusion, we have analyzed strongly correlated 
electron systems on AN. The ground state and thermody- 
namic properties of the Hubbard model have been stud- 
ied using exact diagonalization calculations and quantum 
Monte Carlo. The dependence on the ratio U/t of the 
specific heat, the magnetic susceptibility, the spin-spin 
correlation function, the double occupancy and the one- 
electron transfer support several types of magnetic be- 
havior. We have also studied the magnetic properties of 
the anisotropic spin i Heisenberg model on the AN using 
the real-space RG approach. For ferromagnetic coupling, 
we always observe the existence of ferromagnetism inde- 
pendent on temperature and aniso trop y parameter A. 
As opposed to other structures [Tj, [l5| , the topology of 
the AN favors the ferromagnetic order and the quantum 
fluctuations do not destroy this ordering. For antiferro- 
magnetic coupling, we find a paramagnetic phase for all 
A at finite T. A similar result has been found on AN for 
the Ising model using the transfer matrix method 0. 
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